home *** CD-ROM | disk | FTP | other *** search
/ Shareware Grab Bag / Shareware Grab Bag.iso / 004 / mathstat.arc / M14.BAS (.txt) < prev    next >
Encoding:
GW-BASIC  |  1983-11-13  |  896 b   |  36 lines
  1. 0  CLS:PRINT
  2. 10   PRINT "INTEGRATION: SIMPSON'S RULE"
  3. 30   PRINT "SELECTION: 1=KNOWN FORMULA, 0=UNKNOWN FORMULA";
  4. 40   INPUT S
  5. 49  REM IF FUNCTION IS KNOW ENTER AT LINE 50 DEF FNC(X)='YOUR FUNCTION'
  6. 50   DEF FNC(X)=X
  7. 60   PRINT "LOWER, UPPER LIMIT OF INTEGRAION";
  8. 70   INPUT A,B
  9. 80   PRINT "INCREMENT OF X";
  10. 90   INPUT X1
  11. 95  VR = (B - A) / X1
  12. 100   IF VR <  >  INT (VR) THEN 80
  13. 110   IF S = 1 THEN 150
  14. 120   PRINT "FIRST, LAST VALUE OF F(X)";
  15. 130   INPUT Y1,Y2
  16. 140   GOTO 170
  17. 150  Y1 =  FNC(A)
  18. 160  Y2 =  FNC(B)
  19. 170  C = 0
  20. 180  D = 0
  21. 190   FOR I = 1 TO (B - A) / X1 - 0.5
  22. 200   IF S = 1 THEN 240
  23. 210   PRINT "VALUE OF F(X) AT INTERVAL ";I;" (X=";A + I * X1;")";
  24. 220   INPUT Y
  25. 230   GOTO 250
  26. 240  Y =  FNC(A + I * X1)
  27. 250   IF I / 2 =  INT (I / 2) THEN 280
  28. 260  C = C + Y
  29. 270   GOTO 290
  30. 280  D = D + Y
  31. 290   NEXT I
  32. 300   PRINT "INTEGRAL IS ";X1 / 3 * (Y1 + 4 * C + 2 * D + Y2)
  33. 305  INPUT "Run again (Y/N) ";A$
  34. 307  IF A$="Y" OR A$="y" THEN 10
  35. 310   RUN "MATH"
  36.